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Motivated by the STM experimental data on Bi2Sr2CaCU208+a; which indicate the tunneling conductance 
asymmetry cr(— V) ^ c(V), we report that such a behavior can be explained in terms of the boson fermion model. 
It has been shown in the recent studies, based on various selfconsistent techniques to capture the many-body 
effects, that the low energy spectrum of the boson fermion model is featured by an appearance of the pseudogap 
at T* > T c . We argue that the pseudogap structure has to exhibit a particle- hole asymmetry. This asymmetry 
may eventually depend on the boson concentration. 



1. Introduction 

Recently there was a considerable amount of 
studies of the pseudogap phenomenon observed 
in a variety of experiments on the high temper- 
ature superconductors (HTSC) Q. There are 
two main theoretical interpretations which are 
presently widely considered in the literature: (i) 
the pseudogap as a precursor of the emerging 
pairing fluctuations, and (ii) the pseudogap un- 
derstood in terms of some new (hidden) ordering 
taking place in a vicinity of the superconducting 
phase. Some selective overview can be found e.g. 
in the recent monograph 0] . 

Among theoretical attempts to explain the 
pseudogap effect of HTSC materials there is a 
model of itinerant electrons or holes which coex- 
ist and interact with the local bound pairs (hard- 
core bosons) 0] . It is worth mentioning, that the 
pseudogap has been foreseen within this boson 
fermion (BF) model a long time before the con- 
vincing experimental data became available (see 
the last paragraph in section IV of the Ref. g). 

Pseudogap phase of the BF model is a man- 
ifestation of the pairing-wise correlations which 
start to appear in a system when the transition 
temperature T c is approached from above (the 
precursor type interpretation). On a microscopic 
basis it means that below a certain temperature 
T* fermions start to couple into the pairs. These 



are, however, weakly ordered in phase due to the 
small superfluid stiffness. The phase coherence 
sets in at a sufficiently low temperature T c < T* 
and then the true superconducting transition oc- 
curs. Presence of the incoherent (T* < T < T c ) 
or the coherent (T c > T) fermion pairs is accom- 
panied by either the pseudogap or the true super- 
conducting gap formed around the Fermi energy 
e F . 

In general, some advanced methods of the 
many body theory are required to explore a pseu- 
dogap phase within any model. It is because the 
single-particle and the two-particle correlations 
are then of equal importance. They should be 
properly treated taking account of possible feed- 
back effects between both channels in a controlled 
way. In a context of the BF model such require- 
ments were obtained so far via: a) the selfconsis- 
tent perturbative investigation Pjq], b) the dy- 
namical mean field theory |Q, and c) the flow 
equation study |@|. Some alternative way was 
studied by Micnas et al Q who considered the 
Kosterlitz-Thouless criterion for determination of 
the superconducting transition T c = Tkt and 
identified T* with the mean field estimate Tc ■ 
Authors were able to reproduce qualitatively the 
Uemura type plots T c versus p s . 

In this short paper we extend our previous 
study H to analyze the pseudogap shape and its 
variation with temperature. Direct consequences 
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of the particle-hole asymmetric pseudogap are il- 
lustrated on an example of the STM current con- 
ductance. 

2. Model and the origin of correlations 

We consider the BF model which is described 
by the following Hamiltonian (H 
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where Eq = Ab — 2/j and is the boson energy, 
/i is the chemical potential. The second quantiza- 
tion operators c^J refer to the conduction band 
particles (electrons or holes) and b^' to the com- 
posite hard-core bosons (for instance they can 
represent the trapped electron pairs &j = d^d^). 
Itinerant fermions propagate between the lattice 
sites i and j via the hopping integral tij whereas 
bosons are assumed to be inhnitcly heavy. The 
hard-core property of bosons means that either 
or 1 boson is allowed to occupy a given lat- 
tice site. This local constraint can be formally 
expressed B through the following semi-bosonic 



(l - 2b\bi) 



commutation relations 

and [b h bj} =0 = [&!,&} 

Mechanism of superconductivity and all the 
other forms of correlations in the BF model ([!]) 
are caused by the boson fermion charge exchange 
potential v. By decaying into the fermion pairs, 
bosons gain effectively some mobility. If tempera- 
ture decreases below the critical value T c then (for 
dim > 2) some fraction of bosons gets "frozen" 

into the BE condensate no(T) = ^q = o^q=o^- 

For n (T) ^ 0, fermions are simultaneously 
driven into the broken symmetry superconduct- 
ing state. It can be shown [Q,|8| that the energy 
gap in the superconducting fermion subsystem is 
VyJn {T). 

For temperatures slightly higher than T c there 
exist many bosons which occupy the small mo- 
menta q ~ states. Because of the interaction 



«Ek, q ( & q c k+q/2iC_ k+q / 2T + h.c.) , these boson 
are strongly mixed with the fermion 
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states |n k+q/ 2x) F |n-k+ q /2T) F and thereby the 
life time of fermions might be reduced, especially 
for |k| ~ kp. In consequence, we expect that 
the fermion density of states might be suppressed 
near sp. 

With the on-site boson fermion interaction 
given in (Q) one can generate only the isotropic 
gap/pseudogap. Of course, the HTSC materials 
are characterized by the anisotropic order param- 
eters of the d-wetve symmetry with a possible ad- 
mixture of the s-wave component [jTo| . To cap- 
ture this aspect it is enough to introduce the in- 
tersite coupling Vi t jb\ {ci[Cj^ + CjiCi^ + h.c. when 



both, the superconducting gap |j,|ri|,|12j and the 
pseudogap [|| become anisotropic. Here we only 
discuss the results for the isotropic case but, at 
a price of more difficult numerical computations, 
the same procedure can be easily extended to the 
anisotropic pairing. 



3. The effective spectra 

In order to determine the effective fermion and 
boson spectra of the model (|l]) we utilize the flow 
equation technique proposed by Wegner |l3[ ] . The 
main idea behind is to disentangle the coupled 
boson and fermion subsystems via a sequence of 
canonical transformations H(l) = e~ s ^ He s ^ l \ 
where I is a continuous parameter. We start at 
/ = by putting H BF = H(0), and proceed 
till I = co, when we want to obtain H(oo) = 
H^ ff + Hf ff . All the way, from I = to I = co, 
we adjust the operator S(l) according to the Weg- 
ner's prescription JD|. In practice, disentangling 
of fermion from boson subsystem can be done 
within an accuracy of the order v 3 |8| . To simplify 
the matters we neglect the hard-core constraint 
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and use the pure bosonic relations 
which should be valid for small boson concentra- 
tions tib 

After the disentangling procedure is finished 
we obtain the following structure for the boson 
contribution to the effective Hamiltonian Hfff — 
S q (-^q — ^aO- The initial boson energy Ab is 
thus transformed into the dispersion Eq which 
is characterized by the width proportional to v 2 
and the effective boson mass comparable with the 



mass of fermions ||[§. ^eff P ar t; on the other 
hand, is given as H[ ff = J^w.a ( £ k - m) c L c k<x + 
W Ek, P ,q ^k.p.qC^c^c^Ck+p-qi. Renormaliza- 
tion of £k with respect to the initial dispersion e k 
takes place mainly around k^- There is also in- 
duced the long range fermion-fermion interaction 
C^k, P ,q which has somewhat unusual resonant-like 
character as shown in Figs 7 and 8 of [0 for the 
BCS t/k,-k,q an d for the density-density C/k,q,q 
channels. 

Previously @] we discussed the fermion spec- 
trum only on a basis of the quasiparticle energy 
£k- However, in some cases, a considerable in- 
fluence may also arise from the fermion-fermion 
interactions. These interactions are in princi- 
ple small, |J7k, P ,q| is of the order v 2 , so we can 
treat them perturbatively. The effective disper- 
sion ek a (I) is for T > T c given by 

£k T = £k + jy^^( c W) ■ ( 2 ) 
q 

At T < T c one should also consider the other con- 
tribution coming from the BCS channel E/k.-k.q- 
We restrict our attention only to the normal 
phase (T > T c ). 

In the left h.s. panel of figure |l| we show the den- 
sity of fermion states p{w) = J2k & ( w — £ k)- 
Note, that there appears a pseudogap which 
deepens with a decreasing temperature T. The 
pseudogap structure has a clear particle-hole 
asymmetry at all the temperatures. Asymme- 
try finally disappears at very high temperatures 
T ~ 0.1 (not shown here). The boson den- 
sity of states ^Eq^(^ _ -^q) 1S much less af- 
fected by a varying temperature (see Fig. 4 in 
||). However, upon decreasing T we observe 
(see the right h.s. panel of figure |l|) a consider- 
able redistribution of boson occupancy Nb(&) = 
W E q * (« - #q) fBE (£ q ~ 2/x, T), here f BE is 
the Bose Einstein distribution. By comparing 
both the panels of figure |l| we notice that the 
pseudogap builds up when bosons start populat- 
ing the low energy states £q~o- 

Asymmetry of the pseudogap structure is 
mainly controlled by the boson concentration n E ■ 
Here is a simple argumentation. If the boson en- 
ergy is located in a center of the fermion band 
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Figure 1. The density of fermion states p(oj) and 
the boson occupancy ATg(fi) of the BF model 
with v = 0.1, Ab — and concentrations Uf = 1, 
us = 0.3. Energies u, fl, fi, v and fcsT are all 
expressed in units of the femion bandwidth D. 



(A E — 0), then for the exactly half-filled fermion 
and boson subsystems they both must have sym- 
metric spectra. In particular, the pseudogap 
would then become symmetric too. For the situa- 
tion presented in Fig. |l| we have np — 1, so it can 
only be the boson concentration ub responsible 
for the asymmetry of p(u>). A more detailed anal- 
ysis will be presented in the future publication. 

4. Single particle spectroscopy 

Our results, in particular effects of the particle- 
hole asymmetry, can be well illustrated by cal- 
culating the single particle tunneling current J. 
The differential conductance cr(V) — dJ/dV as 
a function of bias voltage V is a direct probe of 
the density of states below and above the Fermi 
energy. We use the following expression for the 
STM current 

/OO 
djj p(uj) [/jtd(u;,T) 
OO 

-f FD (cj-eV,T)}, (3) 

where Jfd stands for the Fermi Dirac distribu- 
tion. As usually, we neglect the energy w and 
k-dependence of the tunneling matrix |l4| ] . 
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Figure 2. The STM current conductance cr(V) 
as a function of bias voltage V (in units of D/e) 
for the same set of parameters as shown in Fig. 

0- 



In figure || we show the conductance a(V) of 
the STM current (||) obtained for the same set of 
parameters as in Fig. |]. We obtain the negative- 
positive asymmetric characteristics because, at 
low T, the conductance is roughly proportional 
to the density of states p(u>). Our results agree 
very well with the experimental data reported by 
Renner et al Jl5| . Unfortunately, we are unable to 
pass through T c (we solve the flow equations us- 
ing the one dimensional tight binding dispersion 
£k H when Tf m=1 = 0). However, for the realis- 
tic dim > 2, we expect the asymmetry to survive 
even at T < T c as seen experimentally |Q . 

Let us point out the main features of the pseu- 
dogap probed experimentally by the STM con- 
ductance p5| : (i) it is asymmetric, (ii) its magni- 
tude (the peak to peak distance) is almost tem- 
perature independent, and (hi) the pseudogap 
deepens with a decreasing temperature while the 
coherence peaks gradually start to appear. The 
BF model is capable to reproduce all these fea- 
tures (i)-(iii). Some other theoretical concepts 
discussed in the literature to explain a(V), e.g. 
fl4|| and references cited therein, are dealing with 
the physics which microscopically is very close to 
the BF model (Q). 
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